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We calculate the joint distribution P{S, Q) of the scattering matrix S and time-delay matrix 
Q = —ihS'^dS/dE of a chaotic quantum dot coupled by point contacts to metal electrodes. While S 
and Q are statistically independent for ballistic coupling, they become correlated for tunnel coupling. 

We relate the ensemble averages of Q and S and thereby obtain the average density of states at the 
Fermi level. We apply this to a calculation of the effect of a tunnel barrier on the Majorana resonance 
in a topological superconductor. We hnd that the presence of a Majorana bound state is hidden 
in the density of states and in the thermal conductance if even a single scattering channel has unit 
tunnel probability. The electrical conductance remains sensitive to the appearance of a Majorana 
bound state, and we calculate the variation of the average conductance through a topological phase 
transition. 

Contribution for the special issue of Physica E in memory of Markus Biittiker. 


I. INTRODUCTION 

The quantum states of particle and anti-particle exci¬ 
tations in a superconductor (Bogoliubov quasiparticles) 
are related by a unitary transformation, which means 
that they can be represented by a real wave function. In 
this so-called Majorana representation the NxN scatter¬ 
ing matrix S at the Fermi level is real orthogonal rather 
than complex unitary [T]. Since the orthogonal group 
0{N) is doubly connected, this immediately implies a 
twofold distinction of scattering problems in a supercon¬ 
ductor: The subgroup 0+(A^) = SO(A^) of scattering 
matrices with determinant -hi, connected to the unit ma¬ 
trix, is called topologically trivial, while the disconnected 
set 0-{N) of scattering matrices with determinant — 1 
is called topologically nontrivial In mathematical terms, 
the experimental search for Majorana bound states can 
be called a search for systems that have Det S = —1. This 
search has been reviewed, from different perspectives, in 
Refs. IlHSl 

If the scattering is chaotic the scattering matrix be¬ 
comes very sensitive to microscopic details, and it is 
useful to develop a statistical description: Rather than 
studying a particular S, one studies the probability dis¬ 
tribution P{S) in an ensemble of chaotic scatterers. This 
is the framework of random-matrix theory (RMT) [7H9]. 
The ensemble generated by drawing S uniformly from the 
unitary group U(A^), introduced by Dyson in the context 
of nuclear scattering m, is called the circular unitary 
ensemble (CUE). Superconductors need a new ensemble. 
A natural name would have been the circular orthogonal 
ensemble (COE), but since that name is already taken 
for the coset \J{N)/0{N), the alternative name circular 
real ensemble (CRE) is used when S is drawn uniformly 
from 0{N). The RMT of the CRE, and the physical 
applications to Majorana fermions and topological su¬ 
perconductors, have been reviewed recently m- 

The uniformity of the distribution requires ideal cou¬ 


magnetic insulator point contact 



FIG. 1: Andreev billiard on the conducting surface of a three- 
dimensional topological insulator. The billiard consists of a 
conhned region (quantum dot, mean level spacing Jo) with 
superconducting boundaries, connected to metal electrodes 
by a pair of point contacts (supporting a total of A = Ai + 
N 2 propagating modes). A magnetic insulator introduces a 
tunnel barrier in each point contact (transmission probability 
F per mode). A magnetic vortex may introduce a Majorana 
bound state in the quantum dot. 

pling of the scattering channels to the continuum, which 
physically means that the discrete spectrum of a quan¬ 
tum dot is coupled to metal electrodes by ballistic point 
contacts. If the point contact contains a tunnel barrier, 
then P{S) is no longer uniform but biased towards the 
reflection matrix rp of the barrier. The modified distri¬ 
bution PpoissonC^") is known [T2HI6], it goes by the name 
“Poisson kernel” and equals 

-Ppoisson(5') oc Det (1 - (1) 

in the CRE [16] . 

In the present work we apply this result to the scatter¬ 
ing (Andreev reflection) in a superconducting quantum 
dot (Andreev billiard), see Eig. We focus in partic¬ 
ular on the effect of a bound state at the Eermi level 
{E = 0) in the quantum dot, a so-called Majorana zero¬ 
mode or Majorana bound state. In addition to the scat¬ 
tering matrix, which determines the thermal and electri¬ 
cal conductance, we consider also the time-delay matrix 
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Q = —ihS^dS/dE. The eigenvalues of Q are positive 
numbers with the dimension of time, that govern the 
low-frequency dynamics of the system (admittance and 
charge relaxation [MT9]). Moreover, the trace of Q gives 
the density of states and Q and S together determine the 
thermopower [201 III]. 

The joint distribution of S and Q is known for ballistic 
coupling [22142^ . here we generalize that to tunnel cou¬ 
pling. The effect of a tunnel barrier on the time-delay 
matrix has been studied for complex scattering matrices 
[25l|26], but not yet for real matrices. One essential dis¬ 
tinction is that the tunnel barrier has no effect on the 
density of states in the CUE and COE, but it does in the 
CRE. 

The outline of the paper is as follows. The next two 
sections formulate the scattering theory of the Andreev 
billiard and the appropriate random-matrix theory. Our 
key technical result, the joint distribution P(5', Q), is 
given in Sec. |IV| We apply this to the simplest single¬ 
channel case (A" = 1) in Sec. ^ where obtain a remark¬ 
able scaling relation: Eor a high tunnel barrier (trans¬ 
mission probability T <C 1) the distribution P{p\T) of 
the density of states at the Eermi level is described by a 
one-parameter scaling function F{x): 


jF(rp/4) with a Majorana bound state, 

^ |F(4p/r) without a Majorana. 

( 2 ) 

The average density of states in the multi-channel case 
is calculated in Sec. lYg By relating the ensemble averages 
of Q and S we derive the relation 


(p) — (p)ballistic ( 1 ^b] 


(3) 


for a mode-independent tunnel probability T. In the 
CUE and COE the average scattering matrix (S) is just 
equal to tb, so (p) remains equal to its ballistic value 
(p)ballistic: but the CRE is not so constrained. 

Applications to the thermal conductance g and the 


elec trical (Andreev) conductance pA follow in Secs. VII 
and |vml Eor ballistic coupling it is known that P{g) 


is the same with or without the Majorana bound state 

[27] . (This also holds for P{p) [23].) In the presence of 
a tunnel barrier this is no longer the case, but we find 
that the Majorana bound state remains hidden if even 
a single scattering channel has T = 1. The distribution 
of pa: in contrast, is sensitive to the presence or absence 
of the Majorana bound state even for ballistic coupling 

[28] . The way in which P(pa) changes as we tune the 
system through a topological phase transition, at which 
a Majorana bound state emerges, is calculated in Sec. 
El We conclude in Sec.|Xj 

In the main text we focus on the results and applica¬ 
tions. Details of the calculations are moved to the Appen¬ 
dices. These also contain more general results for other 
RMT ensembles, with or without time-reversal and/or 
spin-rotation symmetry. (Both symmetries are broken in 
the CRE.) 


II. SCATTERING FORMULATION 


Fig.[T] shows the scattering geometry, consisting of a 
superconducting quantum dot (Andreev billiard) on the 
surface of a topological insulator, connected to normal 
metal electrodes by point contacts. The Hamiltonian H 
of the quantum dot is related to the energy-dependent 
scattering matrix S{E) by the Mahaux-Weidenmuller 
formula [29] . 

^ l-iwW^E-Hr^W 

^ ’ l+iTrW^E-H)-m ( 4 ) 

= 1 - 2TTiW'‘{E -H + mWW‘')-^W. 

The M X N matrix W couples the M energy levels in 
the quantum dot (mean level spacing (5o) to a total of 
N M propagating modes in the point contact. 

We assume that degeneracies are broken by spin-orbit 
coupling in the topological insulator in combination with 
a magnetic field (perpendicular to the surface). All de¬ 
grees of freedom are therefore counted separately in N 
and M, as well as in (5o. The electron-hole degree of free¬ 
dom is also included in the count, but we leave open the 
possibility of an unpaired Majorana fermion — a coher¬ 
ent superposition of electron and hole quasiparticles that 
does not come with a distinct antiparticle. An odd level 
number M indicates the presence of a Majorana bound 
state in the quantum dot, produced when a magnetic vor¬ 
tex enters m- An odd mode number N signals a propa¬ 
gating Majorana mode in the point contact, allowed by a 
TT-phase difference between the superconducting bound¬ 
aries [3T] . 

The N modes have an energy-independent transmis¬ 
sion probability T^ G [0,1] per mode. If we choose a 
basis such that the coupling matrix W has only nonzero 
elements on the diagonal, it has the explicit form [32] 


fUmn — 

I 12 MSqKyi 

\Wn\ = - 


1 < m < M, 1 < n < A, 


i-r,. 


(5) 


Notice that the tunnel probability T^ determines the 
reflection amplitude G [—1,1] up to a sign. The con¬ 
ventional choice is to take > 0, when Kn = be 

written as 

4 = T(2-r„-2yi-r„). (6) 

J- n 


Alternatively, if < 0 one has K,ri = given by 

K-= T(2-r„ + 2v/rvT;) = iA+ (7) 

J- n 

The two choices are equivalent for ballistic coupling, 
Tn = I = but for a high tunnel barrier En < 1 
one has ^ 0 while ^ oc. The sign change of 
is a topological phase transition [33], which we will ana¬ 
lyze in Section [IX| Eor now we take > 0 for all n, so 

Kn = 
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From the scattering matrix we can obtain transport 
properties, such as the electrical and thermal conduc¬ 
tance, and thermodynamic properties, such as the den¬ 
sity of states. If we restrict ourselves to properties at the 
Fermi level, = 0, we need the matrix S'(O) = S and the 
derivative 


Q 


—ih lim SUE) 

E^O 


dS{E) 

dE 


(8) 


The unitarity of S{E) implies that Q is Hermitian, so 
it has real eigenvalues with the dimension of time. 
The r^’s are called (proper) delay times and Q is called 
the Wigner-Smith time-delay matrix [341136] . The Fermi- 
level density of states p is obtained from Q via the 
Birman-Krein formula [371439] . 


p = lim A In Det S'fE’) = 4_Tr Q. (9) 
2mE^odE ^ ^ 2TTh ^ ^ ' 


For the thermal conductance we partition the modes 
into two sets, N = A^i -h A ^2 5 each set connected to a 
different terminal, and decompose the scattering matrix 
into reflection and transmission subblocks. 


A small temperature difference 5T between the two ter¬ 
minals, at average temperature Tq, drives a heat current 
J = Gthermai<^7^- The thermal conductance Gthermai in 
the low-temperature linear-response limit To,(5T/To ^ 0 
is given by 


5^ — Gthermai/Go — Trtt^, Gq — - ^ ^ • (11) 

The quantum Go is a factor-of-two smaller than in sys¬ 
tems without superconductivity El, due to our separate 
counting of electron and hole degrees of freedom that al¬ 
lows to account for the possibility of propagation via an 
unpaired Majorana mode. 

If we keep the two terminals at the same temperature 
but instead apply a voltage difference, we can drive an 
electrical current. We consider a situation where both 
terminal 2 and the superconductor are grounded, while 
terminal 1 is biased at voltage V. The current I from 
terminal 1 to ground is then given by the Andreev con¬ 
ductance 


h dl 


reer^e + rher'he) 


t _L rr> /r«f 
e( 

= ^ Tr (1 - rr^rV^), 


( 12 ) 


III. RANDOM-MATRIX FORMULATION 


For a statistical description we consider an ensemble 
of quantum dots, each with its own random Hamiltonian 
H. The mean level spacing (5o and coupling matrix W 
are kept fixed. If the wave dynamics in the quantum 
dot is chaotic, the ensemble is fully characterized by the 
presence or absence of certain fundamental symmetries. 
This is the universal framework of random-matrix theory. 

Superconducting systems are characterized by particle- 
hole symmetry, 

H = lF = r,lF*r„ 

^ S'= Ta,5'*r^, Q = r^Q*ra,. 


The Pauli matrices Tx can be removed from the symmetry 
relation by a unitary transformation 


after which we simply have 



H = W = W*, S = S*, Q = Q*. 


(14) 


(15) 


In this so-called Majorana basis the Hamiltonian is real 
antisymmetric, H = iA with A^rn = ^nm — ~^mn- 

If no other symmetries are imposed on the Hamiltonian 
we have the class-D ensemble of random-matrix theory 
mm, with Gaussian probability distribution 

P({^n„^})cx Pexp(-^^). (16) 

The eigenvalues of the antisymmetric M x M matrix H 
come in EE pairs, hence if M is odd there must be a 
nondegenerate eigenvalue = 0 at the Fermi level, in 
the middle of the superconducting gap. This so-called 
Majorana bound state is the hallmark of a topologically 
nontrivial superconductor [42j [43]. If M is even there 
is no level pinned to F^ = 0, and the superconductor is 
called topologically trivial. It is helpful to encode the dis¬ 
tinction in a topological quantum number u that counts 
the number of Majorana bound states, so u equals 0 or 
1 if the superconductor is topologically trivial or non¬ 
trivial, respectively. 

In the scattering matrix the presence of a Majorana 
bound state is signaled by the sign of the determinant. 


Bets'=(-l)^, z/g{0,1}. (17) 

This can be seen directly from the definition 0 in the 
Majorana basis: For M even the matrix H = iA is in¬ 
vertible, so we have 


S = 


IEttW^A-^W 

I-ttIFTA-iIF 


^ Det S = -hi, 


(18) 


in the zero-temperature, zero-voltage limit. In the last 
equality we used the particle-hole symmetry relation t = 
TxEtx at E = 0, where the Pauli matrices act on the 
electron (e) and hole (h) degree of freedom. 


since Det (1 + M) = Det {1—A) if A= —A^. For M odd 
the bound state contributes to the determinant a factor 

_ Det (1 + 

hm ——;-;—^ = —1, 

e^o Det (1 — e~^vv^) 
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for some vector so Bets' = — 1 . 

The class-D ensemble of scattering matrices thus con¬ 
sists of two disjunct sets: The special orthogonal group 
SO(A/') = Oj^{N) of orthogonal matrices with determi¬ 
nant H -1 in the topologically trivial case, and the com¬ 
plement 0-{N) of orthogonal matrices with determinant 
— 1 in the topologically nontrivial case. 


IV. JOINT DISTRIBUTION OF SCATTERING 
MATRIX AND TIME-DELAY MATRIX 

For ballistic coupling (T^ = 1 for all n) the matrices 
S and Q are statistically independent [ 22 ], so they can 
be considered separately. The class-D ballistic scattering 
matrix is uniformly distributed in 0±{N) — uniformity 
being defined with respect to the Haar measure mm- 
This is the Circular Real Ensemble (CRE), the analogue 
for real orthogonal matrices of the Circular Unitary En¬ 
semble (CUE) for complex unitary matrices [SUTO]. 

The class-D ballistic time-delay matrix has probability 
distribution [23] . 

ballistic(Q) oc (DetQ)“^^/^0(Q)exp(-|rHTrQ“^), 

(19) 

where tn = 27rfi/(5o is the Heisenberg time and 0(Q) 
restricts Q to positive definite real symmetric matrices. 
This constraint can be implemented more directly by 
defining 


determines the joint distribution 

P{S, Q) oc (Det (Det 

X exp(-irHTrI]'^(3“^E), (24a) 

^P{S,K) (X (DetI])^exp(-iTrE'^KK'^E). (24b) 

As a check, we can integrate out the time-delay ma¬ 
trix to obtain the marginal distribution of the scattering 
matrix. 



dKP{S,K) oc (DetE)^ 


OK 


-1 


(25) 


The Jacobian evaluates to |33] 


dK 


(DetE)2^-i for K G Rjv,2JV-i, (26) 


and we recover the class-D Poisson kernel^ [16], 


J^Poisson {S) = (Det E) 


l-N 


Det (1 - r^S) 


N-l 


(27) 

The joint distribution (|24|) tells us that S and Q be¬ 


come correlated in the presence of a tunnel barrier. How¬ 
ever, S remains independent of the matrix product 


^0 ^ Q ^ 


(28) 


Q ^ K gRn^2n-i, 


(20) SO that the joint distribution of S and Qq factorizes. 


with ! 


the set of n X m matrices with real elements. 


The distribution (19) is then equivalent^ to the Wishart 
distribution |9| 


-Pwishart(Jf) OC exp{-^ TrK). (21) 


P{S,Qo) — -Ppoisson(^) X ballistic (Qo)- (29) 

The transformation from Q to Qq removes the effect of 
the tunnel barrier on the time-delay matrix (see App. 0 . 


Remarkably, there is no dependence on the topological 
quantum number for ballistic coupling: Q has the same 
distribution irrespective of the presence or absence of a 
Majorana bound state. 

Tunnel coupling is described by a reflection matrix tb 
( from outside to outside) and transmission matrix te 
(from outside to inside). In the Majorana basis these 
are real matrices, parameterized by 

tb = Oi diag(ri,r2,...riv)02, 

tB = Os diag (P}/", ( 22 ) 

0i,02,03 e SO(A^), r„ = 1 e ( 0 , 1 ], 


V. SINGLE-CHANNEL DELAY-TIME 
STATISTICS 


Eor ballistic coupling the distribution (19) implies that 
the eigenvalues 7 ^ = 1 /r^ of Q~^ have the z/-independent 

distibution m 


N 

-Pbaiiistic({7n}) OC exp(-ifH7fe)^(7fe) 

k=l 

e {oa}- 

i<j 


(30) 


As we derive in App. the matrix product 


E = (1 - 


(23) 


The unit step function 0{x) ensures that 7 ^ > 0 for all 
n = 1 , 2 ,...V. 


^ To transform from P(Q) to P{K) multiply by the Jacobians 
\\dQ/dQ-^\\ X \\dKK^/dK\\ = (Detg)^+i x (DetA)^-^ oc 
(Detg)^^/^, for Q = G and K G '^n,2N-i- 


^ In most expressions for the probability distribution we write oc 
to indi cate an unspecified normalization constant. The Poisson 
kernel is normalized, / Ppoisson(S) dS = f dS = 1 with dS 
the Haar measure on 0±(A/’), so we use = instead of oc. 
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FIG. 2: Probability distribution of the inverse delay time 7 
for a single-channel chaotic scatterer, without {u = 0) or with 
(z/ = 1) a Majorana bound state. The histograms are nu¬ 
merical results obtained by generating random Hamiltonians 
(of size M = 40 + z^) with distribution (16). The scattering 
matrix, and hence the delay time, then follows from Eqs. @ 
and The different curves correspond to different trans¬ 
mission probability T of the tunnel barrier. Rescaling with a 
factor , with k dehned in Eq. (32), makes all histograms 
collaps e on to a single curve, in agreement with the analytical 
result ([33I. 


In the single-channel case N = 1 we can use the joint 
distribution to immediately extend this result to ar¬ 
bitrary tunnel probability E = 1 — r^. The scalar S is 
pinned to (— 1 )^, hence 

^ ^ (1 - i-iyrB) ^ I Vk for = 0 , ^ 

yiVTI \l/^ foru = l, ^ ’ 

K = I( 2 -r- 2 \/i-r). (32) 


VI. AVERAGE DENSITY OF STATES 


For ballistic coupling, integration of p = 
(27rfz)“^ 5 ]]^ 7 “^ with distribution ( [M| ) gives the 
average density of states at the Fermi level [23] . 

N 

(p) ballistic = ^^0 ^ _ 2 ^ ^0 = 27rfz/rH, (34) 

for N > 3. The ensemble average diverges for A/" = 1 , 2 . 
To calculate the effect of a tunnel barrier we write 


p= (2^fz)-^Tr(EQoS^), 


(35) 


see Eq. (28), and then use the fact that Qo is independent 
of S and hence independent of E. The average of Qo 
with distribution Pbaiiistic(Qo) is proportional to the unit 
matrix, 


{Qo) = 1 -^(p)ballistic = 1 (36) 

so the average density of states (still for V > 3) is given 
by 

So{p) = ^^Tv{^^^) 

= (e ^ ■ 

(37) 

(This is Eq. ^ from the introduction.) 

It remains to calculate the average of S with the Pois¬ 
son kernel ^7\ . In the Wigner-Dyson symmetry classes 
this average is just rp, but as pointed out in Ref. IT6l this 
no longer holds in the Alt land-Zirnbauer class D. A sim¬ 
ple result for (S) is possible for mode-independent tunnel 
probabilities, F^ = F for all n, see App. 

{S)± = re (l - N-^ ± N-^{1 - ^ (gg) 


[This definition of k corresponds to from Eq. (§.] 
Since k. then appears only as a scale factor, we conclude 
that the single eigenvalue 71 = 7 of Q“^ for = 1 and 
any k G (0,1] has distribution 


where again the + sign corresponds to = 0 (without a 
Majorana bound state) and the — sign to = 1 (with a 
Majorana bound state). 

We thus arrive at the average density of states. 


^ j K^/2exp(-bKfH7) for i/ = 0 , 

^ i/2¥Zh 7 exp(—fK“^fH7) for v = 1 . 

(33) 

The single-parameter scaling P{'y\K^u) = 

zi:i“ 2 ^i^(/^i- 2 ^^^) is tested numerically in Fig. HI by 
drawing random Hamiltonians from the Gaussian class- 
D ensemble ( [T^ . The excellent agreement serves as a 
check on our analytics. 


With p = (27rfz7) ^ the distribution (33) gives the scal¬ 
ing form © of the density of states distribution P(p|r) 
from the introduction, in the tunneling regime F ^ 1 
when ZD = F/4. 


So{p)± 


N 


N-2 



1 ± (1 - r)^/2 


(39) 

plotted in Fig. In the ballistic limit F —>• 1 the depen- 
dence on the Majorana bound state drops out, while in 
the tunneling limit F —> 0 we obtain 


^o{p) = 


2 


N-2 


X 


N-lp 0{T) 
2/F - 1 + 0(F) 


for z/ = 0 , 
for z/ = 1 . 


(40) 


The 1/F divergence of the density of states for z^ = 1 
corresponds to the delta-function contribution from the 
Majorana bound state in the closed system. Without 
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FIG. 3: Ensemble averaged density of states as a function of 
mode-independent transmission probability through the bar¬ 
rier, in the absence {ly = 0) or in the presence = 1) of a 
Majorana bound state. The curves are calculated from the 
analytical expression ( |39[ ). 


the Majorana bound state {u = 0) the density of states 
at the Fermi level remains finite in the F ^ 0 limit, 
but it does remain above the normal-state value oH/Sq. 
This midgap spectral peak is characteristic for a class- 
D superconductor [U [TTJ [42l 031 05] . While in a closed 
system the peak is simply a factor of two, in the weakly 
coupled open system it is a larger factor 2(A^ —1)/(A^ —2), 
which only tends to 2 in the large-A/" limit. The fact 
that the ensemble of open systems does not reduce to an 
ensemble of closed systems in the limit F —0 is due to 
statistical fluctuations that remain important for small 
N. 


VII. THERMAL CONDUCTANCE 


We consider the thermal conductance in the simplest 
case Vi = V 2 = 1 of a quantum dot with single-mode 
point contacts. These are Majorana modes, carrying heat 
but no charge. 

The scattering matrix S G 0±(2) is parameterized by 




cos 0 T sin 0 
sin 0 ± cos 0 


S = 


5'+ if ly = 0, 
S- if u = 1. 


(41) 


The Haar measure equals 

d/i = 7T~^d0^ 0 < < TT, 


(42) 


the same for both 0+ and 0_. 


For the tunnel barrier we take the reflection matrix 
vb = diag(ri,r 2 ), with = \/l ~ > 0. The Poisson 

kernel ^7\ then has the explicit form 


p^^0\ ^ _ 

1 ± rir 2 — (ri ± r 2 ) cos 0 ' 


(43) 


The dimensionless thermal conductance (11) has dis¬ 
tribution 


P±{9) = - [ de S{g - 9)P±{0) 

Jo 

(1 ± rir2)v^rir^ 


= a 


ballistic 


(g) 


(1 ± rir 2)2 - (1 - g){ri ± r 2 Y ’ 


(44) 


where as before, applies to = 0 and P_ to = 1. 
The distribution 

-Pballistic(^) “ / . =; 0 ^ ^ 

is the result m for ballistic coupling (F^ = 1, = 0). 

For identical tunnel barriers, Fi = F 2 = T, this reduces 
to 


-P(^) — -^ballistic (^) ^ 


r( 2 -r) 

r2+4^(i-r) 

1 


if z/ = 0, 
if z/ = 1. 


(46) 


Quite remarkably, the distribution of the thermal con¬ 
ductance for two identical single-mode point contacts is 
unaffected by the presence of a tunnel barrier in the topo¬ 
logically nontrivial case. Fig. is a numerical check of 
this analytical result. 

Notice that the distribution becomes independent 
of z/ if rir 2 = 0. This is a special case of a more general 
result, valid for any ^ 1 ,^ 2 : 


P+{g) = P-ig) if Detre = 0, 


(47) 


in words: The probability distribution of the thermal 
conductance becomes independent of the presence or ab¬ 
sence of a Majorana bound state if the quantum dot is 
coupled ballistically to at least one of the scattering chan¬ 
nels. In other words, ballistic coupling to a propagating 
Majorana mode hides the Majorana bound state. 

The proof is straightforward: If F^o = 1 for one of the 
indices no G {1,2, ...V}, then the Poisson kernel (27) 
is unchanged if we multiply S OiAOjS, with Anm = 
^nm{f — 2(^nno)- [The orthogoual matrix Oi is defined in 
Eq. ([^.] The Haar measure remains unchanged as well, 
and so does the thermal conductance (11). Since Det S = 


—Det (AP), so 0+ is mapped onto 0_, we conclude that 

P+(5) = P-(5). 

This proof for the Poisson kernel extends to the entire 
joint distribution (24) of S and Q: The transformation 


S i-G- OiAOi S has no effect on the matrices Q and E, 
so P(P, Q) remains unchanged. It follows that the prob¬ 
ability distribution of the density of states is the same 
with or without a Majorana bound state if F^ = 1 for at 
least one of the scattering channels. 
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current I is possible for A/^i = 2, when terminal 1 bi¬ 
ased at voltage V has a distinct electron and hole mode. 
We investigate the effect of a tunnel barrier for A/^i = 2, 
N 2 = 1. In the Majorana basis the expression (12) for 
the Andreev conductance reads 


gA = \T^(l-rTyr^Ty), (48) 


with r a 2 X 2 real matrix. 

The scattering matrix S± G 0±(3) can be conveniently 
parameterized using three Euler angles |46], 


_ (R{a) 0\ A 0 \ (R{a') 0\ 
\ 0 ly 1^0 R{0) y V 0 ly 

S- = diag(l,l,-l)5'+, 

where we have defined 

X /cos(a —sin(a 

Ria) = 

^ ^ ysina cos a 


(49) 


(50) 


The Haar measure on 0±(3) is given by [23] 


FIG. 4: Probability distri buti on of the thermal conductance 
g = Gthermai /Go, scc Eq. (11), for a chaotic scatterer having 
two single-mode point contacts with identical tunnel proba¬ 
bilities r. The data points are numerical results for a random 
Hamiltonian (M = 46 -h z^), the curves are the analytical re¬ 
sult (46). In the presence of a Majorana zero-mode {u = 1) 
the distribution is independent of F. 


VIII. ELECTRICAL CONDUCTANCE 

Because a Majorana mode is charge-neutral, no elec¬ 
trical current can be driven for Vi = V 2 = 1. A nonzero 


da = ^^^^dOdada' ^ 


G (0,27r), 6>G(0,7r). (51) 

Because R(a) commutes with the dimensionless 


conductance (48) depends only on the Euler angle 0, 


= 1 ~ cos 0. 


(52) 


In point contact 1 we take a tunnel probability Fi, the 
same for the electron and hole mode, while in point con¬ 
tact 2 we have tunnel probability F 2 for the unpaired 
Majorana mode. Evaluation of the Poisson kernel (27) 
with tb = diag (ri, ri, r 2 ) gives the conductance distri¬ 
bution 


P±{9a) = 


iq/'raKFi - 2 )(±r 2 - 1 ) + (irg - 1 + TM 
[(1 ± r2(5A - 1))2 - (5A - 1 ± r2)2(l - ri)]3/2 


, 0 < flA < 2, r2 = \/l - r2, 


(53) 


plotted in Eig. [^for = 0 (P+), = 1 {P-) and two values of Fi = F 2 = T. 


In the limit r 2 ^ 1, when terminal 2 is decoupled from 
the quantum dot, we recover the result m 


it:;: 


(54) 


independent of Fi. The conductance in this case is 
uniquely determined by the topological quantum num¬ 
ber. 


In the opposite limit r 2 ^ 0 the distribution (53) be¬ 
comes independent of z/. 


P{9a) = 


Fr [2-ri 


(i-ri)5A] 


[i-(i-ri)(flA-W ’ 


if r2 = 0. (55) 


This is a special case of a more general result, for any 


P+(5a) = P-{9a) if DetP2rB = 0, 


(56) 


where V 2 projects onto the modes coupled to terminal 
2. Ballistic coupling, even for a single mode, to terminal 
2 therefore removes the dependence on the topological 
quantum number. 

The proof of Eq. (5^ proceeds along the lines of the 
proof of Eq. (47), with the difference that the transforma¬ 
tion S 1 -^ OiAOjS should leave the upper-left block r of 
S unaffected — otherwise the Andreev conductance (48) 


would change. This also explains why the z/-dependence 
























FIG. 5: Probability distribution of the Andreev conductance 
QA, see Eq. (12), for Ni = 2, iV2 = 1, Fi = r2 = F. The 


data points are numerical results for a random Hamiltonian 
(M = 46 + z^), the curves are the analytical result (53). For 
ballistic coupling (F = 1) the distribution is the same with 
or without a Majorana bound state. In the limit F ^ 0 the 
distribution becomes sharply peaked at = 2z/. 


Majorana bound states 


quantum dot 


metal 

electrode 1 ' 


/' superconductor 

\ 

• • ^ 


ballistic point contact 




nanowire 


FIG. 6: Device to study the topological phase transition in a 
semiconductor nanowire covered by a superconductor. When 
the Zeeman energy of a parallel magnetic field exceeds the 
induced superconducting gap in the nanowire, a pair of Ma¬ 
jorana bound states emerges at the end points. One of these 
is coupled directly to electrode 2, while the other is coupled 
to electrode 1 via a point contact. 


of P(^a) remains for ballistic coupling to terminal 1 [28] . 


IX. MAJORANA PHASE TRANSITION 

The appearance of a Majorana bound state is a topo¬ 
logical phase transition. There is a search for this tran¬ 
sition in a nanowire geometry, see Fig. where it has 
been predicted to occur when the Zeeman energy of a 
magnetic field (parallel to the wire axis) exceeds the gap 
induced by the proximity to a superconductor [48l [49] . 

Because the Majorana bound states emerge pairwise 
at the two ends of the nanowire, the topological quan¬ 
tum number u of the entire structure remains 0 and the 
determinant Det S of the full scattering matrix remains 
+1 through the transition. What changes is the sign of 
the determinant Detr of the reflection submatrix. At 
the topological phase transition Detr = 0, implying a 
perfectly transmitted mode and a quantized peak in the 
thermal conductance [33] . 


We study the effect of the phase transition on the 
statistics of the electrical conductance, measured by con¬ 
tacting one end of the nanowire (terminal 1) to a metal 
electrode at voltage M, while the superconductor and the 
other end of the nanowire (terminal 2) are at ground. 
Terminal 1 is connected to the nanowire via a point con¬ 
tact, thus creating a confined region (quantum dot) with 
chaotic scattering. The minimal dimensionality of the 
scattering matrix S of the quantum dot is 3 x 3: One 
electron and one hole mode connected to terminal 1 and 
one Majorana mode connected to terminal 2. 

To minimize the number of free parameters we assume 
ballistic coupling through the point contact, so the ma¬ 
trix vb in the Poisson kernel (27) is rp = diag(l, l,r 2 ). 
The reflection amplitude r 2 at terminal 2 is tuned 
through zero by some external control parameter typ¬ 
ically magnetic field or gate voltage. Near the transition 
(conveniently shifted to ^ = 0) this dependence can be 
conveniently parameterized by |33] 


r 2 (C) = tanh(CAo)- 


(57) 


(The width of the transition is system dependent.) 
The corresponding coupling constant K 2 in Eq. (§ then 
has an exponential ^-dependence. 


K 2 = I = exp(-2CAo)- 

1 + r2 


(58) 


The probability distribution of the Andreev conduc¬ 
tance (in units of e‘^/h) follows from Eq. (53), 


P{9a) = 5(1 - d)[l + (5 a - l)r2] 


-2 


= i[cosh(^/Co) + (fi-A - l)sinh(^Ao)] 

0 < SA < 2. 


-2 


(59) 


The delta-function limits (54) are reached for ^ ^ Too 


(keeping z/ = 0, because the entire system is topologically 
trivial). Right at the transition, at ^ = 0, the distribution 
is uniform in the interval 0 < < 2. The average 

conductance varies through the transition as 


(5a) = 1 


+ 


e/Co 


tanh(^Ao) sinh^(^Ao) 


(60) 


Fig.0 


X. CONCLUSION 

In conclusion, we have investigated a variety of observ¬ 
able consequences of the fact that the scattering matrix 
of Majorana fermions is real orthogonal rather than com¬ 
plex unitary. Of particular interest is the identification 
of observables that can detect the sign of the determi¬ 
nant, since Bets' = — 1 signifies the presence of a Ma¬ 
jorana bound state. The obvious signal of such a zero¬ 
mode, a midgap peak in the density of states [l2j|43], is 
broadened by tunnel coupling to the continuum. We find 
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FIG. 7: Variation of the ensemble-averaged Andreev conduc¬ 
tance in the geometry of Fig. as the nanowire is driven 
through a topological phase transition (controlled by a pa¬ 
rameter which can be thought of as the deviation of the 
magnetic held from the critical held strength). The curve is 
the result (60) for a single-channel ballistic point contact. 


that the peak remains hidden if the coupling is ballistic 
(unit transmission) in even a single scattering channel. 
The thermal conductance is likewise insensitive to the 
presence or absence of a Majorana bound state, but the 
electrical conductance retains this sensitivity when the 
coupling is ballistic. 

These results for the effect of a tunnel barrier on the 
midgap spectral peak are derived from the distribution 
P(5', Q) of scattering matrix and time-delay matrix un¬ 
der the assumption of chaotic scattering, due to disor¬ 
der or due to irregularly shaped boundaries. The ap¬ 
propriate ensemble in the absence of time-reversal and 
spin-rotation symmetry has symmetry class D in the 
Alt land-Zirnbauer classification [1]. Chiral symmetry 
would change this to class BDI, in which multiple zero- 
modes can overlap without splitting [32]. The effect 
of chiral symmetry on the joint distribution P{S,Q) is 
known for ballistic coupling [24] — but not yet for tun¬ 
nel coupling. This seems a worthwhile project for future 
research. 
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Appendix A: Joint distribution of scattering matrix 
and time-delay matrix 

We calculate the joint distribution P{S, Q) of scat¬ 
tering matrix and time-delay matrix in the presence of 
a tunnel barrier, starting from the known distribution 
P{So,Qo) without a barrier [22l[23|. The application in 
the main text concerns symmetry class D, but for the 
sake of generality and for later reference we give results 



Altland-Zirnbauer 

Wigner-Dyson 


D 

Dill 

C 

Cl 

A 

AI 

AH 

a 

-1 

-1 

2 

1 

0 

0 

0 

p 

1 

2 

4 

2 

2 

1 

4 

to/rn 

1 

1 

1 

2 

1 

2 

1 

1 

1 

degeneracy d 

1 

2 

2 

2 

1 

1 

2 


TABLE I: Parameters that appear in the distribution of the 
scattering matrix and time-delay matrix, for each of the 
Alt land-Zirnbauer and Wigner-Dyson symmetry classes. No¬ 
tice that a different set of indices o', govern the energy level 
statistics m- The degeneracy factor d refers to the Kramers 
degeneracy of the scattering channels and the delay times, ig¬ 
noring uncoupled spin bands. (The energy levels may have a 
different degeneracy.) 


for all four Altland-Zirnbauer [T] symmetry classes D, 
Dill, C, Cl, as well as for the three Wigner-Dyson laisnj 
symmetry classes A, AI, AIL The symmetry indices that 
distinguish the ensembles are listed in Table [l] see Ref. 
ITT] for an overview of this classification. 

The unitary matrix S and the Hermitian matrix Q are 
real in class D, complex in class A, and quaternion in 
class C. We consider these three symmetry classes with¬ 
out time-reversal symmetry first, and then include the 
constraints of time-reversal symmetry in classes Dill, Cl, 
AI, AII.^ 


1. Broken time-reversal symmetry 

Without the barrier So is independent of Qo and uni¬ 
formly distributed, 

P{So, Qo)dfi{So)dfi{Qo) = P{Qo)dfi{So)dfi{Qo). (Al) 

The differential d/i indicates the Haar measure for the 
unitary matrix So and the Euclidean measure for the 
Hermitian matrix Qo- The ballistic time-delay matrix 
distribution is given by [22| |23] 

P(Qo-i) K (Det'Qo-i)“+^^/2 

X e(Qo)exp(-i^ioTr'Qoi), (A2a) 

4^ P{Qo) oc 

X e(Qo) exp(-i^io Tr' (A2b) 

Degenerate eigenvalues of Qo are counted only once in 
Tr' and Det'. In terms of the degeneracy factor d from 
Table HI this can be written as 

Det'Qo = (DetQo)i/t IV'Qo = I IVQo. (A3) 


^ The seven symmetry classes in Table^do not exhaust the tenfold 
way classification of random-matrix theory: There are three more 
chiral classes m (labeled AIII, BDI, CII) that require separate 
consideration m- 
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D 

C 

A 

AI 

Cl 

All 

Dill 


broken time-reversal symmetry 

preserved time-reversal symmetry 


6S = S'^dS = -5S^ 

<55 = UUSU^ = -<55+ = <55 



III 

U = U° = ayU'^ay 

SSnm 
(n / m) 

go 

13 = 1 

qoCFo +iq ‘ cr 
13 = 4. 

a + ib 
/5 = 2 

ia 

(3=1 

iaax + ibaz 
P = 2 

iqoao + q ■ cr 
(3 = 4 

aax + bcfz 
13 = 2 

SSnn 

0 

0 + 1 = 0 

iq • a 

0 + 1 = 3 

ib 

0 + 1 = 1 

ia 

0+1 = 1 

iaax + ibTz 

0 + 1 = 2 

iqoao 

0 + 1 = 1 

0 

0 + 1 = 0 


Q = -ihSUS/dE = Qt 

Q = -im\dS/dE)U^ = Ql =Q 

Qnm 

^0 

qoao -\-iq CT 

a-\- ib 

a 

acTo + ibay 

qoao -\-iq cT 

acTo + ibay 

Qnn 

^0 

qocro 

a 

a 

aao 

qocro 

aao 


TABLE II: Characterization of the scattering matrix differential SS and of the time-delay matrix Q. All coefficients Qn, ci,b are 
real, and ai is a Pauli matrix. The symmetry indices /3 and a1 count, respectively, the number of degrees of freedom of the 
off-diagonal and diagonal components of the anti-Hermitian matrix SS. The off-diagonal elements of the Hermitian matrix Q 
have /3 degrees of freedom, while the diagonal elements have one single degree of freedom in each symmetry class. 


The channel number N also does not include degenera¬ 
cies, so the total number of eigenvalues of Qo is d x N. 
The characteristic time to differs from the Heisenberg 
time tn by a numerical coefficient,^ see Table [ij 

Insertion of the barrier, with unitary scattering matrix 


determine the Jacobian [44] 


Js 


d/i(5') _ ^SSo^"^ 

dfi{So) ~ SSo 

(Det'EE’^)^^“b/3/2+i+a^ 


(A8) 


Sb = 



(A4) Eq. (A6) also implies the relation between the time- 

delay matrices. 


transforms So into 


Q = SQoEI' ^dQ = SdQoSi' + 0{dS). (A9) 


S — Tb + ^B'^o) 

4^ 5o = - 5 Vb)(1 - 

Variations of S and So are related by [9] 


The off-diagonal elements of the Hermitian matrix Q 
have P degrees of freedom, the diagonal elements have 
one single degree of freedom in each symmetry class. The 
Jacobian is then given by 


SUS = T.{SldSo)T.\ S = (1 - (A6) 

The differentials 


Jq = 


dfjj(^Qo^ 


'EdQo'^^ 


= (Det'EE’^) 


dQo 

(Ar-l)/3/2+l 


(AlO) 


SS = SUS, SSo = SldSo, (A7) 

are anti-Hermitian matrices, = —SS. The number 
of degrees of freedom of the off-diagonal elements are 
given by p and the number of degrees of freedom of the 
diagonal elements by 1 + a. As summarized in Table [TI| 
real matrices (class D) have a = — 1, /3 = 1, complex 
matrices (class A) have a = 0, /3 = 2, and quaternion 
matrices (class C) have a = 2, /3 = 4. These parameters 


The joint probability distribution P(5', Q) now follows 
upon division of P{So, Qo) by the product of Jacobians, 


P(P,Q) 


P{Qo) 

JsJq 

P(go)(Det' 


(All) 


Substituting Qo = E ^ 

rive at the joint distribution 


into Eq. (A2) we thus ar- 


^ We define tn = 27r/i/(5o, with (5o the mean spacing of nondegener¬ 
ate levels. The ratio to/tn then equals the degeneracy of energy 
levels divided by the degeneracy of delay times m- It is unity 
in all symmetry classes except C and Cl, where the delay times 
have a Kramers degeneracy that the energy levels lack m- 


P{S,Q) OC (Det'SSl')'^^/2(j)g^/Q)-3/3iV/2+/3-2-a 

X e(Q)exp(-i/?ioTr'EtQ-iE). (A12) 

The class-D result p4| ) from the main text follows for 
a = —1, /3 = 1, to = th, d = 1, = E"’’. 
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2. Preserved time-reversal symmetry 

Time-reversal symmetry equates the scattering matrix 
to its transpose in class AI and Cl and to its dual 

in class All and DHL (The dual of a matrix is = 
(TyS^cTy.) We use a unified notation S = where the 
tilde indicates the transpose or the dual, whichever is 
appropriate for that symmetry class. The symmetry S = 
S allows for the “square root” factorization 

S = uu = S, (A13) 


with unitary U. 

The time-delay matrix in these symmetry classes is 
constructed such that it satisfies the same symmetry, 


Q = —ih lim = Q. (A14) 

E^o dE 

This redefinition of Q differs from Eq. by a unitary 
transformation, so the delay times are not affected. 

The ballistic Qo and So are agai n ind epen dent [2^ |23] , 
distributed according to Eqs. (Al) and (|A2| ) with the ap¬ 
propriate values of a and P from Table [Ij These numbers 
now count the diagonal and off-diagonal degrees of free¬ 
dom of the symmetrized differential 


5S = EdSU\ 


(A15) 


constrained by 5S = —5S\ 5S = 5S. The matrix ele- 
ments of 6S are imaginary in class AI (a = 0, /d = 1), i 
times a quaternion^ in class All (a = 0, /3 = 4), and of 
the form aaxEbaz with a, b imaginary in class Cl (a = 1, 
/3 = 2) and a, b real in class Dill {a = —1, P = 2). 

The elements of the Hermitian matrix Q are real 
in class AI, quaternion in class All, and of the form 
acjo + ibay with a, b real in both classes Cl and DHL 
The off-diagonal elements of Q have the same number of 
/3 degrees of freedom as (55', but the diagonal elements 
have only a single degree of freedom irrespective of a. 
All of this is summarized in Table HD 

The symmetrization of the differential modifies the re¬ 
lation ( |A6| ), which now reads 

6S = U^U^SSoUo'E^U^, S = (1 - (A16) 

The relation ( |A9[ ) between Q and Qo is similarly modified 
by the symmetrization, 

Q = C/SC/o^Qot/oStf/t. (A17) 


Beca use t he m atrices I/, Uq are unitary, the Jacobians 
(A8) and (AlO) are unchanged. 


Js = (Det' 5]5]t)(Af-i)/3/2+i+a^ (AIS) 

Jq = (Det' 5]5]t)(iv-i)/3/2+i_ (Aig) 


We thus obtain the joint distribution 

P{S,Q) OC (Det'SSl')'^^/2(-j)g^/Q)-3/3iV/2+/3-2-a 

X e(Q)exp(-i/3toTr'l]1'C/1'Q-i?7E). (A20) 


3. Poisson kernel 


The marginal distributio n of the sc atter ing matrix re¬ 
sulting from the Jacobians (A8) and (A18) is 


P{S) = J dQP{S,Q) = l/Js 

= (Det' 5]5]t)-(^-i)/3/2-i-« 


Det'(l - rps) 
|Det'(l -4S')|2 


(JV-l)/3/2+l+« 


(A21) 


including the normalization constant. This formula com¬ 
bines the known expressions for the Poisson kernel^ in 
the Wigner-Dyson ensembles m and in the Altland- 
Zirnbauer ensembles m- 


The present analysis confirms that Eq. (A21) holds 
without modification in the two symmetry classes D and 
Dill that support Majorana zero-modes, depending on 
the sign of the determinant Det 5 = ±1 in class D and 
the sign of the Pfafhan Pf {icfyS) = ±1 in class Dill [52] . 
As a check, we can take N = when 5 = ±1 in class D 
and 5 = Tao in class DHL The ± sign determines the 
presence or absence of a Majorana bound state (twofold 
degenerate in class Dill). Since there is only a single el¬ 
ement in the ensem ble w e should have P{S) = 1, which 
is indeed what Eq. dA^ gives for = 1, O' = —1. 


Appendix B: Calculation of the ensemble-averaged 
scattering matrix 

1. Symmetry class D 


According to Eq. (37), the effect of a tunnel barrier 


on the average density of states follows directly once we 
know the average scattering matrix. Simple expressions 
can be obtained if we assume that the tunnel probabilities 
are mode-independent, T^ = T for n = 1, 2 ,... A. 

The scattering matrix of the barrier has the polar de¬ 
composition 


5b = 




0 O3 



® The name “Poisson kernel” applies strictly speaking only to 
the Wigner-Dyson ensembles, when vb = f SP{S)dS. In the 
Altland-Zirnbauer ensembles the average scattering matrix dif¬ 
fers from re, see Ref. 1161 and App. 


^ A quaternion has the form aoao + iaiax + ia 2 (Jy + iasaz, with 
four real coefficiients an ■ The matrix cro is the 2x2 unit matrix. 
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with Oi, 02 , 03,04 G SO(A^) . T he block structure cor¬ 
responds to Eqs. ( [^ an d (|A4| ), in particular, vb = 
VT^T^Oi02 . From Eq. ( |A5f we obtain the scattering 
matrix S of the quantum dot, 


S = Oi 


VT^ + ru{i + VT^u)-^ O 2 , (B2) 


2. Symmetry class C 

In a similar way we can derive the average scattering 
matrix in other symmetry classes. We give the result for 
class C. The matrix U then varies over the unitary sym- 
plectic group Sp(2A/'). Ref. [53l gives the required average: 


in terms of a matrix U = O 45 'o 03 that is uniformly dis¬ 
tributed in 0±{N). 

Because the average of for any power p = 1, 2 ,... 
is proportional to the identity matrix, we may write the 
average of S in the form of a power series, 

{S)± = (^1 - ^ E(-ir(l - Tr/^^{TrUP)±^ . 

(B3) 

If the average of U would be over the entire unitary 
group, then all terms in the power series would vanish 
and we would simply have {S) = tb- But averages over 
orthogonal matrices do not vanish, in the nontrivial way 
calculated^ by Rains [53] : 


{TiU^)c 


— 1 if p < 2N and even, 
0 otherwise. 


(B6) 


Eq. (B3) still holds with the factor 1/N replaced by 
1/2A^. We thus find the average scattering matrix in 
class C, 


(S)c=rB(l+^[l-{l-my (B7) 

For N = 1 this gives {S)c = (1 + r/2)rB, in agreement 
with Ref. Iini 

The average density of states in class C follows from 
Eq. where we again account for the doubling of the 
dimensionality N 2N: 


(TiUn ^ + if p- TV = 0,2,4,... 

^ 2 10 otherwise. 

(B4) 


We substitute Eq. (B4) into Eq. (B3) and sum the geo¬ 


metric series, to arrive at the average scattering matrix 


{p)c — (p) ballistic ( 1 [('^) C ^b] 


N 


{N + I)(5o 


l-i^[l-(l-rn). (B8) 






(B5) 


used in Sec. [V^ to obtain the average density of states in 
class D. 


In the second equation we have substituted the ballistic 
class-C result from Ref. [23l The tunneling limit T ^ 0 
gives a vanishing density of states, 

NT 


consistent with the class-C result for a closed system [42l 

143] . 

^ For the record, we note that Eq. ( |B4| ) differs from the formula in 
Ref. 1531 bv a minus sign (zb instead of =f)- 
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